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The inverse problem with Lemaitre-Tolman-Bondi (LTB) universe models is discussed. 
The LTB solution for the Einstein equations describes the spherically symmetric dust-filled 
spacetime. The LTB solution has two physical functional degrees of freedom of the radial 
coordinate. The inverse problem is constructing an LTB model requiring that the LTB 
model be consistent with selected important observational data. In this paper, we assume 
that the observer is at the center and consider the distance-redshift relation Da and the 
redshift-space mass density fi as the selected important observational data. We give Da and 
\x as functions of the redshift z. Then, we explicitly show that, for general functional forms of 
Da(z) and fJ,(z), the regular solution does not necessarily exist in the whole redshift domain. 
We also show that the condition for the existence of the regular solution is satisfied by the 
distance-redshift relation and the redshift-space mass density in /1CDM models. Deriving 
regular differential equations for the inverse problem with the distance-redshift relation and 
the redshift-space mass density in ACDM models, we numerically solve them for the case 
(Qmo, &ao) = (0.3,0.7). A set of analytic fitting functions for the resultant LTB universe 
model is given. How to solve the inverse problem with the simultaneous big-bang and a 
given function Da(z) for the distance-redshift relation is provided in the Appendix. 

§1. Introduction 

Anti-Copernican universe models are widely discussed in recent years as alter- 
natives to the standard homogeneous cosmology with dark energy components. One 
of the simplest ways to construct an anti-Copernican model is to solve the inverse 
problem, which is the reconstruction of the universe model from observational data. 
Although the isotropy of the universe around us has been confirmed with high ac- 
curacy by the observation of the cosmic microwave background (CMB), this does 
not automatically imply homogeneity of the universe. Thus, in solving the inverse 
problem, we may assume that the universe is spherically symmetric around us. In 
addition, we usually assume that the universe is dominated by cold dark matter, 
that is, by dust. The spherically symmetric dust-filled spacetime is described by the 
Lemaitre-Tolman-Bondi (LTB) solutionP^2} The LTB solution has three arbitrary 
functions of the radial coordinate approximately corresponding to the density profile, 
the spatial curvature and the big-bang time perturbation, with one of them being a 
gauge degree of freedom representing the choice of the radial coordinate. These arbi- 
trary functions may be determined by requiring that the resulting LTB universe be 
consistent with selected important observational data (e.g., the distance-redshift re- 
lation). However, we should note that it is not apparent at all if these three functions 
have sufficient degrees of freedom to fit all of the important observational data. 

In the inverse problem, two approaches have been mainly considered. One is 
proposed by Mustapha, Hellaby and Ellis^ in 1998. In this paper, we refer to this 
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work and the approach used there as MHE and the MHE approach, respectively. In 
MHE, the angular diameter distance D\ and the redshift-space mass density \x are 
given as functions of the redshift z to fix two physical functional degrees of freedom 
in LTB universe models. 

Recently, many authors have succeeded in solving the inverse problem by the 
MHE approachP^SJ anc i constructed the LTB universe model whose distance-redshift 
relation and the redshift-space mass density agree with those in the concordance 
ylCDM model. The other approach is proposed by Iguchi, Nakamura and Nakao in 
2002. We refer to this work as INN in this paper. In the INN approach, one of 
the conditions is given by the distance-redshift relation as in the MHE approach, 
and simultaneous big-bang or a uniform curvature function is chosen as the other 
condition. However, they could not go beyond z ~ 1.6 owing to a technical prob- 
lem. In 2008, Yoo, Kai and Nakao succeeded in constructing the LTB model whose 
distance-redshift relation agrees with that of the concordance ylCDM model in the 
whole redshift domain and which has uniform big-bang timeP^J 1 Independently from 
MHE and INN, Celerier solved the inverse problem analytically at small redshifts 
z <C 1 in the form of the Maclaurin series in 1999. 1 ™ There are also several works 
on the inverse problem {EHED 

In the inverse problem, one of the difficulties happens at the point with the 
maximum angular diameter distance. At this point, differential equations for the 
inverse problem become apparently singular. Related discussions can be found in 
Refs. B|)-f7|) JT6|) . The point with the maximum angular diameter distance is a regular 
singular point of the differential equations in the INN approach. In Ref. 10), this 
singularity has been resolved using a shooting method to solve the differential equa- 
tions for the inverse problem. On the other hand, there is no conclusive illustration 
on how to resolve this apparent singularity in the MHE approach. The main purpose 
in this paper is to explicitly show how this singularity can or cannot be resolved in 
the MHE approach. 

In this paper, we solve the inverse problem by the MHE approach using a differ- 
ent formulation from MHE. We find that, for general functional forms of D\(z) and 
fi{z), the regular solution does not necessarily exist in the whole redshift domain as 
pointed out in Ref. U]). Then, we represent the necessary and sufficient condition 
for the existence of the regular solution in terms of Dj±(z) and fj,(z)F^ We also show 
that this condition is satisfied by the distance-redshift relation and the redshift-space 
mass density in TlCDM models. Deriving regular differential equations for the in- 
verse problem with the distance-redshift relation and redshift-space mass density in 
ylCDM models, we numerically solve them for the case (I?mo> &ao) = (0.3,0.7) and 
compare our result with those in previous works J^~^ We propose a set of analytic 
fitting functions for the resultant LTB universe model. We also explain how to solve 
the inverse problem by the INN approach with simultaneous big-bang in Appendix 
ICl We use the unit system given by c = G = Hq = 1 throughout this paper, 
where c, G and Hq are the speed of light, Newton's constant and Hubble constant, 
respectively. 
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§2. Equations 



2.1. LTB dust universe 

As mentioned in the introduction, we consider a spherically symmetric inhomo- 
geneous universe filled with dust. This universe is described by an exact solution of 
the Einstein equations, known as the LTB solution. The metric of the LTB solution 
is given by 



ds z 



2 , (d r R(t,r)) 2 2 2 



-dr z + R\t,r)dQ\ 



(24) 



1 — k(r)r 2 

where A;(r) is an arbitrary function of the radial coordinate r. The matter is dust 
whose stress-energy tensor is given by 

T» v = pu^u u , (2-2) 

where p = p(t, r) is the mass density, and u a is the four-velocity of the fluid element. 
The coordinate system in Eq. (|2-ip is chosen in such a way that = (1, 0, 0, 0). 
The area radius R(t, r) satisfies one of the Einstein equations, 



{d t Rf 



-k(r)r 2 + 



2M(r) 



-k(r)r + 



2 m(r)r 3 



R w 3R ' 

where M(r) is an arbitrary function related to the mass density p by 

<9 r M _ \r 2 m{r) + ^d r m{r) 



47T/3 



and we have defined 



R 2 d r R 



M(r) 



R 2 d r R 

m(r)r 3 
6 ' 



Following Ref. [T7|) . we write the solution of Eq. (|2-3p in the form, 

R(t,r) = {§M{r)) x l z {t - t B {r)f' z S{x) = rm{r) l ' z {t - t B {r)) 2 / 3 S{x) 



x = k(r)r 2 



2 ( t - t B (r) 
6M(r) 



2/3 



k(r)m(r)- 2 / 3 (t-t B (r)) 2 / 3 , 



(2-3) 

(24) 

(2-5) 

(2-6) 
(2-7) 



where t B (r) is an arbitrary function that determines the big-bang time, and S(x) is 
a function defined implicitly as 



S(x) 



COSh yJ — T] — 1 



-(sinh yj—f] 



-v 



)2/3 



6 1 / 3 (sinh yf—r\ — y/—n) 2 / 3 ' 
1 — cos ^/fj 

6 1 / 3 (y / r/ — sin y^) 2 / 3 ' 



6 2/3 

{y/rj-sin ^T]) 2 / 3 



6 2/3 



for x < . 



for x > . 



and 5(0) = (3/4) 1 / 3 . The function S(x) is analytic for x < (tt/3) 2 ^ 3 . Some charac- 
teristics of the function S(x) are given in Refs. fTU|) and [T7|) . 
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2.2. Basic equations 

As shown in the preceding subsection, the LTB solution has three arbitrary 
functions: k(r), m(r) and t#(r). O ne of them is a gauge degree of freedom for 
rescaling of the radial coordinate r. We fix the remaining two functional degrees of 
freedom imposing the following physical conditions. 

• The angular diameter distance is given by Da(z) as a function of the redshift. 

• The redshift-space mass density is given by /j,(z) as a function of the redshift. 
To determine t-s(r), k(r) and m(r) from the above conditions, we consider a 

past directed outgoing radial null geodesic that emanates from the observer at the 
center. This null geodesic is expressed in the form 

t = t(z), (2-9) 
r = r(z). (2-10) 

We assume that the observer is located at the symmetry center r = and observes 
the light ray at t = to. To fix the gauge freedom to rescale the radial coordinate r, 
we adopt the light-cone gauge condition such that the relation 

t = t -r (2-11) 

is satisfied along the observed light ray. 

Then, basic equations to determine ts, k and m are given as follows: 

1. Null geodesic equations 

The null geodesic equations in the LTB solution are given as 

. . dt o r R , . 

(1 + ^ = -^7r (242) 

{L + Z) dz~ d t d r R ■ [2L6) 

One of these equations can be derived from the other one and the null condition. 
Therefore, one of them is sufficient under the null condition. 

2. Null condition 

By virtue of the light-cone gauge condition, the null condition on the observed 
light ray takes the very simple form of 

d r R =y/l- kr 2 . (2-14) 

3. Distance-redshift relation 

As mentioned, we assume that the angular diameter distance is given by D^(z). 
Then, we have 

R = D A (z). (2-15) 

4. Redshift-space mass density 

We define the total mass of a shell between z and z + dz as 



(2-16) 
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Then, we consider the redshift-space mass density defined by 



M dr _^ t _ 

v = ir-R2 f i = pit- 2-17 

A-kDI dz 



As mentioned, we assume that the redshift-space mass density is given as a 
function of the redshift. Thus, we have 

dr , . , 
p— = n{z). (2-18) 

2.3. Rewriting the basic equations as differential equations 

We can regard r(z), m(z), k(z) and £3(2) as mutually independent functions of 
the redshift z. Basic equations to determine these functions are (j2- 13|) . (|2-14|) . ([2-15P 
and (jMS]). 

Combining Eqs. (|2-4p and (|2-18p . we find 

dr dm / : — -Dl 

3m— +r— = 24vr^Vl - kr 2 —^ = A, (249) 

dz dz r 1 



where we have used the relation 



Differentiating Eq. (|2-6p with z, we have 



„dr dm dk dr dD\ 

B— + C— + D— + E— = —±, 2-21 

dz dz dz dz dz 



where 



B = m x / 3 (i - t B ) 2/3 S - ^rm^it - t B y 1/3 (S + xS'), (2-22) 

C = l rm - 2/3 (t - t B ) 2/3 (^ - 2xS'), (2-23) 

L> = rm~ 1/3 (t-tB) 4/3 S', (2-24) 

£ = -^rm^it - t B )~ 1/3 (S + xS'). (2-25) 

From Eqs. ()2-6|) and (|2-14p . we have 

F^- 0^-2^-^ = 0, (2-26) 

dz az dz 



where 



F= VT^kr 2 -m 1 ^{t-tvf l3 S. (2-27) 
For dtd r R, we find the expression 

d t d r R = Gd r m + Hd r k + Id r t B + J, (2-28) 



6 



Yoo 



where 



G = d t C = ^rm- 2 '\t - t^y 1 / 3 ^ (2-29) 

H = d t D = \rm- l/z {t - t B ) 1/3 (2S' + xS"), (2-30) 

I = d t E= irmV3(t - t B )- 4/3 ^ , (2-31) 

J = ~m l / 3 (t-t B )- 1/3 {S + xS'). (2-32) 
3 

Then, Eq. (|2- 13|) can be rewritten as 

jdr | c dm | R dk | ^dfe _ Vl-kr 2 
dz dz dz dz 1 + z 

Solving Eqs. $2lM . (f2^2TD . (12^261) and (12^33]) for dr/dz, dm/dz, dk/dz and 
dr/dz, we have 



dr _ 1 dD A 
dz ~ B + F dz ' 

G??n 1 / . „ dr\ If. 3m dD&\ 



(2-34) 

■ A-Sm^ )=-[A- ) , (2-35) 

cterV dzj r\ B + F dz J ' y 1 



where C = C/r, D = D/r, E = E/r, G = G/r, H = H/r and I = I/r, and we have 
used the relation 

DI-EH=-. (2-38) 

§3. Solution near the center 

Before numerically solving Eqs. (|2-34|) - (|2-37|) . we need to identify the boundary 
conditions for the differential equations at the center. Since the center is the regular 
singular point of the differential equations, the boundary conditions must be appro- 
priately chosen for a regular solution. For this purpose, we consider the Maclaurin 
series expansion of the functions near the center as follows: 

r = r lZ + ^r 2 z 2 + 0{z 3 ), (3-1) 

m = m + m\z + 0(z 2 ), (3-2) 
k = k + hz + O(z 2 ), (3-3) 
*B = t B1 z + 0{z 2 ). (3-4) 
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Hereafter, each value at (t,r) = (to,0) is denoted by the subscript 0. 
As shown in Ref. [TOlh imposing the regularity at the center, we find 

R = r + 0(r 2 ), (3-5) 

m = mo + O(r), (3-6) 

m Q = 87rpo := 8irp(t , 0), (3-7) 

, 2/3,-2/3 i , m /q n 

k = x m ' t = -1 + — , (3-8) 

So = mo V V /3 - (3-9) 

Since S satisfies 



S + x S' = ^±-3x for x <{^ 2 '\ (3-10) 



S + xS' = ^^-3x for x>(£j 2/ \ (3-11) 



we find 



S + x S' = -m 1/3 tl /3 , (3-12) 



where we have assumed 



^mo = ^ < 1 < x < < (^) 2/3 . (3-13) 

We assume the following expansion form of the angular diameter distance and 
the redshift-space mass density: 

D A = z+^D A2 z 2 + 0(z 3 ), (3-14) 

8TTfj, = do + SiZ + Oiz 2 ). (3-15) 

Then, from zero-th order of Eq. (|2-34p . we can find 

ri = 1. (3-16) 

Equation ()2-35j) can be expanded as 

dm mn — <5n „„ , s ,„ _ _ s 

— = -3—^ - + 0(1). 3-17) 

dz z 

Therefore, we have to impose the condition 

m = 5 Q (3-18) 
so that dm/dz is finite at the center. Assuming the above relation, we can find 

r2 = l + D A2 (3-19) 
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from the second order in Eq. (|2-34p . In the same manner as in the above, we obtain 

m 1 = -^((2 + D A2 ) 5 - 5 1 ), (3-20) 

fci = 2 + D A2 - 1(2 + Da2)<5o + j, (3-21) 
-2(3to + 2)dg + 3{5 1 + 12)t 6 Q - 3D A2 ((5 - 6)i + 2)5 + 18tfi(t„ - 1) 

from the 1st orders of Eqs. (|2-35p - (|2-37p . where we have used the equation 

x(2SS" + S' 2 ) + 5SS' + - = (3-23) 

to reduce the order of differentiation of S. Then, we can use these expressions near 
the center instead of solving Eqs. ()2-34[) — (|2-3T|) . If D A (z) and fj,(z) are analytic at 
the center, we can derive higher order expressions for Eqs. (|3-1|) - (|3-4|) from Eqs. 

It should be noted that all the initial valuables at the center are determined 
by the input valuables So and Hq = 1 from the regularity. We do not have any 
additional degree of freedom to put a boundary condition for the differential equa- 
tions differently from the situation in Ref. 10). We also note that the normalization 
Hq = 1 implicitly determines the value of to through Eq. (|3-9j) . 

§4. Regularity at the point with the maximum distance 

To obtain a physically reasonable solution, the right-hand side of Eq. (|2-34j) 
must be positive definite. Since the sign of dD A /dz changes at the point with the 
maximum angular diameter distance, at which dD A /dz = 0, B + F must be at this 
point and the sign of B + F must change. As shown in Appendix |XJ we can derive 
the equation 

B + F = J ±-(l-Q(z)), (4-1) 

where 

Q(z) = 4tt / (1 + z)n{z)D A {z)dz. (4-2) 
Jo 

Therefore, the input functions D A {z) and fj,(z) must satisfy 

Q(z m ) = 4vr / (1 + z)n(z)D A (z)dz = 1 (4-3) 
J o 

for the regularity at z = z m , where dD A /dz\ z=Zm = 0. This condition has been 
pointed out in RefJD). 

Let us consider the angular diameter distance and the redshift-space mass den- 
sity in a homogeneous and isotropic universe model with ]?mo , ^ko and J?aq . In this 
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case, the angular diameter distance is given by 
D A (z) = D A CDm(z) : = 



lh Jo H^ d * for ^K Q = 0, 

V-^KO fo H^ dz \ for tt K0 ^0, 

(4-4) 



1 sin 



where 



Hacbm = V(l + z) 3 Q M o + (1 + z) 2 n KQ + Q A o- (4-5) 
The redshift-space mass density is given by 



/ \ / \ 3 (1 + z) 3 (dD AC DU , D AC t>m\ 
= WCDM z := t- Ouo—^= -j + — 



ylCDM 

Then, the following equations are satisfied: 



(4-6) 



d_ 

dz 



2 tt dD ACBM 



(1 + zYHacbm 



3(1 + 2 



2 



,3 



dz 

jV 1 + f? K0 (l + *) 2 ^ C dm = L 1 + ^ VkoDacdm for ^K0 + 0, (4-8) 



ylCDM 



■1?M0-C/1CDM, (4-7) 



[(1 + 2) 2 £» 2 1CDM] = AdCDM for ^K0 = 0. (4-9) 

Using these equations, we can find 

Q(z) =4tt (1 + z)nAC-DM(z)D A CT)M(z)dz 
Jo 

= [ ^M0(1 + zfQ^^l + % (1 + ^) 2 ^icDM 

= jT -| [(1 + ,) 2 ^ CDM ^f^] Cfe = 1 - (1 + ,) 2 i/,CDM^f^lO) 

for J?ko 7^ 0. We can also derive the same equation for the case 1?ko = 0> as 
pointed out in Ref. [8]). Therefore, the condition (|4-3p is automatically satisfied for 
homogeneous and isotropic universes with i?M0> ^ko and i?Ao- in this case, we have 

B + F = (l + z)H ACDM ^^. (4-11) 



Then, Eq. (|2-34p can be replaced by 

dr 



dz (l + z)H 



ylCDM 



(4-12) 



The set of differential equations (|4-12|) . (|2-35p - (|2-37p does not have any singular 
point other than the center. 

On the other hand, if we consider a dark energy component with the equation 
of state p = wp (w / -1) instead of the cosmological constant A, the condition 
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w 

Fig. 1. Q(Zm) as a function of w. Da and fi are given by those in homogeneous and isotropic 
universe models with Qmo = 0.3 and i?xo = 0.7. 



cannot be satisfied. We numerically calculated the value of Q(z m ), with ,J?mo = 0.3 
and fixo = 0.7, where J?xo is the density parameter of the dark energy component. 
The result is shown in Fig. [TJ This result shows that there is no regular solution for 
the inverse problem with and fi for a homogeneous universe model with a dark 
energy component with w ^ 1 in the whole redshift domain. In other words, the 
existence of the regular solution for the inverse problem cannot be guaranteed for 
general input functions D\(z) and p,{z) without the condition (|4-3j) as pointed out 
in Ref. HJ. 



§5. Numerical results 



5.1. Numerical integration 

In this subsection, we consider the distance-redshift relation Dj± and the redshift- 
space mass density fx for a homogeneous and isotropic universe with (4?mo> ^xo) = 
(0.3, 0.7), where J?xo is the density parameter for the dark energy component whose 
equation of state is given by p = wp. The distance-redshift relation satisfies the 
following second-order differential equation and boundary conditions JIH1>>Q11> 



dz 2 



d2DA - {Th + dJ P) d ^t- t + a + «>fto(i + *r- 2 ) , 

(5-1) 

^a(0) = dD ^ 



dz 



1, (5-2) 



H = v / (l + z)3^ M0 + (l + z)«'+i^ X0 . (5-3) 



where 



The redshift-space mass density is given by the same expression as Eq. 

We numerically integrated Eqs. (|2-34j) — ()2-37|> . Equation (|4-12|) is used instead 
of Eq. (|2-34p for the case w = —1. Near the center, solutions can be given in forms 
of the Maclaurin series, as shown in Appendix |Bl We use these expressions near the 
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center instead of solving the differential equations. As shown in Fig. [2j we can find 
a regular solution in the whole redshift domain only for w = —1. For w = —1, r(z), 
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Fig. 2. Radial coordinate along the past directed null geodesic as a function of redshift for each 
value of w. 



m(z), k(z) and £3(2) are depicted in Fig. [3J m, k and ts are depicted as functions 
of r in Fig. SI The energy density p on the t = to time slice is depicted as a function 
of circumferential radius R(to,r) in Fig. [SJ This hump type profile has already been 
found in Ref. /£) and consistent results can be seen in Refs. [8]) and [9]). Our result 
for w = — 1 is also consistent with the result in these previous papers. 

5.2. Gauge transformation and fitting functions 

We propose fitting functions for the given solution. Before that, for convenience, 
we perform a gauge transformation r = r{r). Hereafter, the "~" represents the 
quantities in the new gauge. We consider the new gauge defined by 

~ 4-7T n 

M(r) = M(r) = —r 3 p . (5-4) 
3 
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Fig. 3. r, k, m and tB as functions of z. 
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Fig. 4. k, m and ts as functions of r. 
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Fig. 5. p(to,r) as a function of R(to,r) for the numerically obtained model and the model with 
analytic fitting functions. 



Then, the relation between r and r is given by 

~3 = m(r) r3 
m 

The curvature function k(r) and the big-bang time is (j) are given by 

^(r)? 2 = k(r)r 2 , 
t~B(7) = t B (r). 

We propose analytic fitting functions for k and £b as 
7 



(5-5) 



(5-6) 
(5-7) 
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+ 0.728893F - 0.634917F 2 + 0.303959F 3 - 0.073768f 4 



0.00878216, 
r + 0.3039 



*Bflt(r) = -0.21319F - 0.013605? 2 + 0.000925931F 3 + 



0.298252r 
r + 1.5439 



(5-9) 
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Fig. 6. fc(f), fcflt(r), tsir) and tBfit(r). 

The differences in the angular diameter distance and redshift-space mass density 
between the ylCDM model and LTB model with k(r) = k^ t {r) and t-^ir) = t-Q^ir) 
are less than 1%, as shown in Fig. [7J 
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Fig. 7. Differences in the anguiar diameter distance and redshift-space mass density between the 
ylCDM model and LTB model with k(r) = fc fit (r) and t B (r) = t Bfit (r). 
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§6. Summary 

As part of the summary, we describe a theorem and two corollaries based on the 
analysis in the previous sections after giving three definitions for convenience. 

Definition 1 

We say that -Da (z) is reasonable as a function for the distance-redshift relation 
if it satisfies the following: 

• Dx{z) is a C 1 function on z > and is piecewise smooth on z > 0. 

• D A (z) > for z > and Z? A (0) = 0. 

• ^L =i /tfo=i. 

• ^ = 0atz = %. 

• D A (z) is C 2 and ^fr" < at z = z m 
. ^ ^ for z ^ z m . 

Definition 2 

We say that //(z) is reasonable as a function for the redshift-space mass density 
if it satisfies the following: 

• fi(z) is finite and positive definite. 

• /i(z) is a piecewise smooth function on z > 0. 

Definition 3 

We say that an LTB universe model is observationally regular if the LTB universe 
model is regular on the past light-cone of the observer at the center except for the 
big-bang initial singularity. 

Theorem 



For a set of a reasonable angular diameter distance D\(z) and a reasonable 
redshift-space mass density fi(z), there exists the observationally regular LTB uni- 
verse model whose distance-redshift relation and redshift-space mass density for the 
observer at the center agree with D A {z) and ^(z), respectively, if and only if Da(z) 
and fi(z) satisfy 



4vr / (l + z)fj,(z)D A (z)dz = l. (6-1) 
Jo 

Proof 
Since 

B + F = ^-L- (l - 4vr J (1 + z)n(z)D A (z)dzj , 

Eq. ([EH) is equivalent to B + F = 0. We can find that (1 + z){B + F) is a 
monotonically decreasing function of z, and B + F can vanish only once in z > 0. 
If£? + i ? /0atz = z m , since dD A / dz\ z=z = 0, we have dr/dz\ z=z = from Eq. 
(|2-34p . Then, from Eq. (|2-18p . we find that p must be infinite at z = z m for a finite 
value of ^(z m ). Therefore, B + F must vanish at z = z m for the existence of the 
regular solution. That is, Eq. (|6T|) is a necessary condition for the existence of the 
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regular solution. If the condition (|6Tp is satisfied, from L'Hopital's rule, we have 

d 2 D A /dz 2 



dr 
dz 



> 0. 



Then, we can find an observationally regular LTB solution solving Eqs. (|2 34|) (|2 37|) . 
Q.E.D. 

Corollary- 1 

For an observationally regular LTB universe model in which fi(z) and D A {z) 
are reasonable, the distance-redshift relation D A {z) and fi(z) satisfy 4tt / Q Zm (l + 
z)/J,(z)D A (z)dz = 1. 

Proof 

Since Eqs. (|2-34p - (|2-37p are applicable for any LTB universe model, 4ir J* Zm (l + 
z)[x{z)D A (z)dz = 1 must be satisfied so that an LTB universe model is observation- 
ally regular. Q.E.D. 

Corollary-2 

There exist observationally regular LTB universe models whose distance-redshift 
relation and redshift-space mass density coincide with those in /1CDM models. 
Proof 

As shown in the text, -D^icdm^) and hacdm(z) satisfy the condition ()6Tp . Then, 
this corollary is an immediate consequence from the above theorem. Q.E.D. 

The condition (|6-ip has been derived in Ref. H|) at the first time. Without 
this condition, we cannot obtain the LTB universe model beyond z = z m . We 
demonstrated this using D\(z) and fi(z) for flat homogeneous and isotropic universe 
models with the dark energy component whose equation of state is given by p = wp. 
Ifw 7^ —1, we cannot obtain the regular LTB solution beyond z = z m by solving the 
differential equations $M$-i$Z2Q. 

In £j5j we have obtained the LTB universe model that realizes the distance- 
redshift relation and the redshift-space mass density in the TlCDM model with 
(]?MO)^Ao) = (0.3,0.7). Then, we introduced analytic fitting functions for this 
numerically obtained LTB model. The LTB model with these analytic fitting func- 
tions realize the distance-redshift relation and the redshift-space mass density in the 
ylCDM model within 1% accuracy. 
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Derivation of Eq. 



Before the derivation, we provide some useful equations. 
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Prom the expression for B + F, we can find 

B + F = —d t R + d r R. (A-l) 
Differentiating Eq. ()2-3p with t, we have 



mr 



3 



afl? = - (A-2) 

* 6F 2 



Differentiating Eq. ()2-14[) with z along the light-cone, we have 
d „ _ d 



— d r F = — y/l-kr 2 ^ d 2 R = d t d r R + d r \/l - kr 2 , (A-3) 

where we have used the relation 

d dr dt _ dr .„ _ . ,. 



Differentiating Eq. (|2-3p with r and dividing both sides by ^/T^hr 2 , we have 



mr 3 d t d r Rd t R d r (mr 



6F 2 VI - kr 2 6FV1 - kr 2 
First, differentiating Eq. (|A-lj) with the redshift z, we have 

iL(fl + F ) = ^ (d*F + t 2 i? - 2d t d r F) . (A-6) 



Using Eqs. (|A-2p and (|A-3|) . we can rewrite the above equation as 

, P> 4- Fl = — liLs^ - l:r 2 - — ' 

dz 

Then, using Eqs. (|A-5[> and (|2-14p . we have 



>^) = ^(^VT^-^-^). (A-7) 



d - (B + F) = d 4- ( mR (d t R - d r R) dr{mr3 



dz y dz VVl - kr 2 6RVl - kr 2 

dr ( d t d r R {B + F) _ d r (mr 3 ) \ ^ 



dz V y/l - kr 2 6FV1 - kr 2 

Using Eqs. J2H3I), an d (EiMD, we can find 

±{B + F) = -^±^-A^D K 
dz 1 + z 

& — ((1 + z)(F + F)) = -4tt(1 + z)iiD A . (A-9) 
Since (F + F)| 2=0 = 1, we have 

B + F = ^-j— (l - 4vr y (1 + z)fiD A dz^j . (A-10) 
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Appendix B 

Series Expansion of the Solution 

We consider the distance-redshift relation and the redshift-space mass den- 
sity n for homogeneous and isotropic universes with (J?mo> ^xo> ^ko = 1 — ^mo — ^xo) 
as input functions for the inverse problem given by a set of equations, Eqs. (|2-34|) 
- (|2-37|) . For the solution near the center, we write functions r(z), m(z), k(z) and 
ts(z) in the Maclaurin series as follows: 

r = nz + K 2 z 2 + K^z 3 + ^r A z 4 + 0(z 5 ), (B-l) 

m = tuq + m\z H — m 2 z 2 H — ra^z 3 + 0(z 4 ), (B-2) 
2 6 

k = k + hz + \k 2 z 2 + -k 3 z 3 + 0{z 4 ), (B-3) 
2 o 

t B = taiz + \t B2 z 2 + \t m z 3 + 0(z 4 ). (B-4) 

2 6 

Substituting these expressions into Eqs. (|2-34]) - (|2-37p . we can find the values of each 
coefficient as shown below. For notational simplicity, we use y := 3w + 1. 

r(z) = z + j(-f2 U o - ^xoy - 4)z 2 

+ 24 ( 3n M0 + (6^xoy + 8)^ M0 + ^xo(y((3^xo - 2)y + 8) - 4) + 24) z 3 

+ y^(-15^mo " 9(5^xoy + 4)^o " 3(I?xo(y(3(5^xo - 2)y + 26) - 8) + 24)^ M0 
+ f2 X0 [y ((3(6 - 5^xo)^xo - 4)y 2 - 42fl X0 y + 28y + 24tt X o - 92) + 56) - 192) z 4 
+ O (z 5 ) , (B-5) 

m(z) = 3f2 MO + + |-^mo(-50^mo + 24^ X0 + Wf2 xo y + 60)z 2 

2 oU 

/ 2 2 2 2 

+ ™^mo(25^ mo - 12^xo ^mo + 20i7xo2/ ^mo - 50I2 M o - 5^ X o2/ + 5 ^xo2/ + 22i7 X o 

oU 

- 12^ X0 y + 20^ X0 y + 10) z 3 + O (z 4 ) , (B-6) 



k(z) = -1 + J2 M0 + ( ^M0 - ^7T^ - 1 ) Z 



+ (-100/^0 + 2(^ X0 (25y + 36) + 80)tt M0 - 5^xo(y((8 - 3^xo)y + 16) + 16) - 60) z 2 



-±-(l00^ + (^xo(50y - 72) - 200)^ 



1 

4~80 

+ (^xo(y(50y - ^xo(65y + 72) + 190) + 148) + 100) Q M0 

- 5^xo(y(y((3(^xo " 2 )^xo + 4)y + 4) + 24) + 8))z 3 + O (z 4 ) , (B-7) 
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tn(z) 



(I?xo(3t - 2)y)z 




MO 



+ ((92 - 81i )^xo + 32(6t - 5).f2xo - 48t + 32) ^ M0 
+ 2(12t + 3^xo (-26*o + ^xo(28t - 23) + 20) - 8))y 3 

- 10(i?Mo — 1)(— 24i?xo^Mo + 27i?xo*o^Mo ~~ 18*o^mo + 12/2 MO + 30i?xo^Mo 

- 54^ X o*o^MO - 30t ^M0 + 20I2 M o + 84I2 X o - 108^x0*0 + 48i - 32)y 2 

+ 4(tt M0 - 1)(-45*o^mo + 30^ - 36I?xoio^Mo + 120* ^mo " ^mo ~ 40^ X o^mo 
+ 126l?xoio^Mo - 255t ^Mo + 90^mo + 60^ X0 - 120^xo*o + 180t - 80)y 

- 8(I?mo - I) 2 (l8t ^ M0 " Sto^MO + 18^ M0 - 48t + 8))z 3 + O {z 4 ) . (B-8) 



Here, we show how to solve the inverse problem using the INN approach with 
^b(^) = for a given angular diameter distance D\(z). This has been carried out 
in Ref. [T0|) by solving a set of four differential equations parametrized by an affine 
parameter on the null geodesic. Here, we do not use the affine parameter but the 
redshift z as the independent variable. In this procedure, the number of differential 
equations is reduced to three. 

C.l. Basic equations 

Basic equations are derived by replacing the condition (|2- 18[) with ts(z) = 0. 
Then, we have the following three coupled first-order differential equations: 



Appendix C 

Inverse Problem with t-g{r) = 



dr 

dz 



1 dD A 

B + F dz 



(C-l) 




(C-2) 



(C-3) 
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C.2. Solutions near the center 

We perform the Maclaurin series expansion as (j3- lf l— f)3- 14|) . Substituting these 
expressions into Eqs. (|C-1[) - (|C-3[1 . we find 

r(z) = z + X -(D k2 + l)z 2 + ^(4D A2 + 2D A3 + m )z 3 + O (z 4 ) , (C-4) 
m(z) = m + r "° (7 "° + ^ 2 Q ( +" 6) 9 / °L" 6 18t ° + 9> + O (z 2 ) , (C-5) 

V ' 3 3((mo + 6)t - 6) 

Then, we find that the parameter mo can be freely chosen differently from the MHE 
approach^ We use this deg ree of freedom to guarantee B + F — Oat z — z m . 

C.3. Numerical procedure to solve the differential equations 

For illustrative purposes, we rewrite Eqs. (|C-ip - (|C-3p in abstract forms as 

dr 1 dD A 



dz X(z,r,m,k) dz 



(C-7) 



fz = z { z ^ kd iz d i^)- (c ' 9) 

Then, the procedure to solve Eqs. (|CTp - (|C-3p is summarized in follows: 

1. We determine a trial value for mo- 

2. We solve the differential equations from the center, where, near the center, we 
use the expressions (|C-4p - (|C-6p instead of solving the differential equations. 

3. We stop integrating the equations at z = z^ < z m and read off the values of r, 
m and k at this point. We label these values as r^,-, mb- an d kb_. 

4. We determine trial values for m m and k m , where the subscript "m" denotes the 
value at z = z m . 

5. We numerically solve the equation X(z m , r m , m m , k m ) = for r m . 

6. At z = z m , Eq. ()C-7p can be rewritten as 



dr 

dz 



d 2 D A 



^-{z, r, m, k, dr/dz, dm/dz, dk/dz) 



(C40) 



Then, Eqs. (|C-8p - (|CT0p can be regarded as simultaneous equations for dr/dz \ 
dm/dz\ z=Zm and dk/dz\ z=Zm . We numerically solve them. 
7. We integrate the differential equations (|C-ip - (jC-3p from z = z m backward. 



*' We accept a nonvanishing first derivative of the density at the centerP^'^ 
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8. We stop integrating the equations at z = zy, and read off the values of r, m and 
k at this point. We label these values as r b+ , m b+ and fc b+ . 

9. We define deviations 5r, 5m and 5k as follows: 



5r 
5m 
6k 



r b _-r b+ , (C-ll) 
m b _-m b+ , (C-12) 
h--k h+ . (C-13) 



10. Operations 1-9 give the deviations 5r, 5m and 5k for given values of mo, m m 
and k m . That is, we can define 5r = 5r(mo,m m ,k m ), 5m = 5m(mo,m m ,k m ) 
and 5k = <5/c(mo, m m , fc m ) as functions of mo, m m and k m using operations 1-9. 
Then, repeating operations 1-9 and using the Newton-Raphson method, we 
numerically solve the equations 

5r(m ,m m , k m ) = 0, (C-14) 
5m(m Q , m m , k m ) = 0, (C-15) 
5k(m ,m m , k m ) = (C46) 

for mo, m m and k m . Eventually, we have functional forms of r(z), m(z) and 
k(z) in z < z m without discontinuity at z = z b . 

11. Finally, we solve the differential equations (|CTp - (|C-3p from z = z m forward 
with initial conditions given above. 

C.4. Results 

We performed operations 1-11 in the previous subsection using the distance- 
redshift relation in the ylCDM model with (fiyio, ^Ao) = (0.3,0.7). The result is 
consistent with that in Ref. [10]) . We performed the same gauge transformation as 
that given in Then, we find a fitting function for k(r) as^ 



0.545745 2.22881 r , _ 

feet = ; — ■ (C-17 

0.211472 + V0.026176 + f ( .807782 + ^0.026176 + r) 2 



The angular diameter distance in the LTB universe model with k{r) = k^r) is 
depicted in the Fig. [5J The deviation of the distance from that in the ylCDM model 
is within 0.1%. 
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